Land subsidence analysis along high-speed railway based on EEMD-Prophet method

Environmental changes and ground subsidence along railway lines are serious concerns during high-speed railway operations. It is worth noting that AutoRegressive Integrated Moving Average (ARMA), Long Short-Term Memory (LSTM), and other prediction methods may present limitations when applied to predict InSAR time series results. To address this issue, this study proposes a prediction method that decomposes the nonlinear settlement time series of feature points obtained through InSAR technology using Ensemble Empirical Mode Decomposition (EEMD). Subsequently, multiple Intrinsic Mode Functions (IMFs) are generated, and each IMF is individually predicted using the Prophet forecasting model. Finally, we employ an equal-weight superimposition method to combine the results, resulting in the prediction of the InSAR settlement time series. The predicted values of each component are subsequently weighted equally and combined to derive the final prediction outcome. This paper selects InSAR monitoring data along a high-speed railway in inland China and uses the proposed method and ARMA and Prophet models to carry out comparative experiments. The experimental results show that compared with the ARMA and Prophet models, the method in this paper improves the root mean square error by 58.01% and 32.3%, and increases the mean absolute error by 62.69% and 33.78%, respectively. The predicted settlement values generated by our method exhibit better agreement with the actual InSAR monitoring values.

the HLSTM method that accounts for spatial heterogeneity, thereby mitigating the influence of spatial heterogeneity on the prediction results of large-scale land subsidence.However, using the above method combined with InSAR to predict the subsidence field needs to establish a complex model and consider the influence of various factors.If some factors are not taken into consideration, it may cause significant errors in the prediction results and other shortcomings 12 .
The Prophet model 13 is a time series forecasting model based on an additive model.In 2017, Taylor et al. proposed a new signal-processing method and released an open-source software package.This model employs a generalized additive model to conduct Bayesian curve fitting of nonlinear trends on time series data, thus offering remarkable flexibility and robustness 14 .The EEMD 15 is based on the empirical mode decomposition 16 (EMD) method, which effectively addresses the modal aliasing of the IMF obtained via EMD decomposition by superimposing Gaussian white noise 17 .The core concept is to perform multiple decompositions on different scales of the input data.By averaging the outcomes of each decomposition, the randomness and problems with local extreme values in EMD decomposition are reduced, and the decomposition results are more stable.
Aiming to address the aforementioned issues, this study proposes a prediction method that decomposes InSAR nonlinear time series features before predicting them.Specifically, the settlement time series obtained by InSAR is firstly decomposed into intrinsic mode functions using the EEMD method, which can characterize local features at different times.Then, the Prophet model is used to predict each intrinsic mode function.Finally, the predicted results of each intrinsic mode function are combined with equal weights to enhance the prediction accuracy.

Method Main idea
The method proposed in this paper is based on the SBAS-InSAR technology to predict the settlement time series characteristics of InSAR settlement.SBAS-InSAR was originally proposed by Berardino et al. 18 .This technology optimizes the combination of images acquired at different times, and its main process is shown in the literature 19 .
In our approach, we first perform EEMD decomposition of the SBAS-InSAR time series results.Subsequently, we make separate predictions for each of the decomposition results, and finally, we apply a novel method involving equal-weight superposition.Our method combines the strong adaptability, fast prediction speed, and high prediction effect of the Prophet model with the EEMD method's characteristics of not requiring prior conditions.Figure 1 illustrates the flow chart of our method for predicting InSAR nonlinear subsidence time-series data, and the main steps are as follows: (1) Conduct SBAS -InSAR processing on the study area to obtain the InSAR subsidence time series results.

SBAS-InSAR Processing Workflow
SBAS-InSAR is a monitoring technology based on traditional D-InSAR technology, which is further deepened to obtain the time series characteristics of surface deformation.The algorithm selects N + 1 SAR images, sets an appropriate spatiotemporal baseline consistent with PS-InSAR, generates N small baseline set differential interference pairs, and performs phase unwrapping after generating differential interference patterns.During this process, we set Multilook Number to 4 where Range Looks is 4 and Azimuth Looks is 1.This reduces the impact of coherence spots and achieves a good balance between signal-to-noise ratio and spatial resolution.In addition, in order to eliminate atmospheric noise, Goldstein filtering is used, and its window size is set to 5*5.The selection of Unwrapping Method Type is Minimum Cost Flow, where Unwrapping Decomposition Level is 1 and Unwrapping Coherence Threshold is 3.These parameter choices were guided by careful consideration of the trade-offs between computational efficiency and preserving deformation signals in the presence of atmospheric effects.
It is assumed that within the ordered time (t 1 , t 2 , ...t n ) of the interest area, N + 1 SAR images of the same orbit are acquired, and the auxiliary images are registered to the common main image according to the interference conditions, and then the auxiliary images are registered to the common main image according to the interference conditions.The complex phase solution is performed on the image, and M differential interference patterns are finally generated.Among them, M needs to meet the following conditions: During the interference process, in the jth differential interference pattern generated by the two SAR images acquired at t 1 and t 2 ( t 1 > t 2 ) as the main and auxiliary images, the interference phase of any pixel can be expressed as: In the formula, is the wavelength, d(t 1 , x, r) and d(t 2 , x, r) are the cumulative deformation quantities, �φ j top is due to the reference DEM Terrain phase error caused by inaccuracy, �φ j atm is the error phase caused by atmospheric delay, �φ noise is the error phase caused by the influence of noise.After removing the terrain phase error, the error caused by atmospheric delay, and the error caused by the influence of noise, through simplified formula (3), the interference phase can be expressed as: To obtain the sedimentation sequence, the average phase velocity can be expressed as: In the formula, the phase value of the jth interference pattern can be expressed as: The deformation of the time series calculated by the SBAS-InSAR algorithm can be expressed as: In the formula, A is the MXN matrix, and 80 is the vector representing the interference phase value.The SBAS-InSAR technical data processing flow chart is shown in Fig. 2.

Prophet prediction model based on EEMD
The Prophet prediction method, based on EEMD, combines the advantages of both techniques.It can better capture the nonlinear, non-constant, and periodic changing characteristics of the time series while effectively removing noise and interference through decomposition.This approach enhances prediction accuracy while preserving the important information of the original signal.In line with the proposed EEMD-Prophet prediction method, the original signal of the InSAR monitoring point is firstly decomposed into several IMF and residual components using the EEMD method, specifically: (1) Add a set of Gaussian white noise n(t) to the InSAR settlement time series signal feature Q(t) to obtain an overall signal X(t) (2) Decompose X(t) by EMD first to obtain the corresponding IMF of each order, and its formula is expressed as: (3) Add a set of different white noise signals n j (t) to the InSAR settlement time series feature Q(t) to obtain another overall signal X j (t) .After performing EMD decomposition on X j (t) , obtain the corresponding IMF components of each order using its formula expression: (1) www.nature.com/scientificreports/(4) Repeat the above steps, and calculate the mean value of each IMF obtained as the IMF of the final signal, and its formula is expressed as: where IMF i (t) denotes the IMF for each specific order, whereas res i (t) corresponds to the associated residual component.N signifies the total count of iterations involving white noise addition.t is the time for obtaining InSAR results, and i, j are the number of IMFs obtained by decomposition.
According to the principle of the Prophet model, each IMF obtained by decomposing the characteristics of the InSAR subsidence timing signal by the EEMD method is adaptively decomposed into: In Eq. ( 11), the Prophet model decomposes each IMF signal time series feature y imf (t) into individual ele- ments: the time series data trend item g imf (t) , the oscillation period item s imf (t) , the holiday mutation term h imf (t) , and the residual term ε t .The decomposed trend item g imf (t) represents the changing function of the nonlinear growth component within the time series.When analyzing land subsidence time series, the trend item is commonly captured using a logistic regression function 20 .which can be represented as: In Eq. ( 12), k represents the settlement growth rate, l denotes the displacement, and m represents the upper limit of the settlement trend.As time t increases, g imf (t) gradually converges towards m.However, the time series of surface subsidence is subject to changes, requiring consideration of the dynamic subsidence growth www.nature.com/scientificreports/rate k and the adaptive offset.Additionally, the changing growth rate needs to be taken into account.Hence, Eq. ( 12) is revised as: This Eq. ( 13) describes the variation of the trend component g imf (t) in the ground subsidence time series with respect to time t.In this context, the growth rate k and the offset l are dynamically changing with time, influenced by the parameters α(t), δ, and γ.This dynamic modeling allows the trend to adapt over time to better fit the observed data.The symbol t stands for time, signifying the time point.The symbol α(t) is a function related to time t.This function represents parameters that vary over time.The symbol δ is a parameter vector, signifying the variation in the growth rate k over time.The symbol γ is a parameter vector, representing the variation in the displacement l over time.
The oscillating period item s imf (t) includes seasonal periodic oscillations and small trend changes, and the Fourier series construction of the settlement time series data can obtain s imf (t) , the structural form can be expressed as: In Eq. ( 14), T represents the length of the timing cycle, while 2n denotes the anticipated number of cycles in the model.Typically, the land subsidence time series does not exhibit a holiday mutation term h imf (t) , and thus, its influence is disregarded.The residual term ε t adheres to a normal distribution and can be interpreted as unexpected random noise and varying trends.The technical approach of the EEMD-based Prophet prediction method is illustrated in Fig. 3.

Study areas and data
The Taiyuan-Jiaozuo high-speed railway is situated in the inland region of eastern mainland China.It serves as a vital transportation link connecting Jiaozuo City in Henan Province to Taiyuan City in Shanxi Province.This railway line commenced operations on December 12, 2020, spanning a total length of 358.764 km and encompassing 13 intermediate stations.Designed for a maximum speed of 250 km per hour, the route traverses a diverse terrain comprising plains, hills, mountains, rivers, and other intricate landscapes, each presenting unique geological characteristics 21 .Notably, the railway passes through certain coal mine goaf areas, introducing additional uncertainties to train operation safety.Given the stringent requirements for track stability in high-speed rail systems, even minor surface deformations pose substantial risks to the overall operation.Figure 4 illustrates the geographical location of the research area, specifically focusing on the section of the study section spanning from Jinzhong Station to Taigu East Station along the Taijiao high-speed railway.
The data utilized in this study comprises ESA Sentinel-1A data, obtained during the satellite image period spanning from January 09, 2021, to October 07, 2022, encompassing a total of 40 scene data.The corresponding timestamps for each SAR image are presented in Table 1.The dataset is based on C-band observations with an ascending orbit direction and an incident angle of 39.27°.The imaging mode employed is the interferometric wide mode IW, with VV polarization.The pixel resolution is 5m in the range direction and 20m in the azimuth direction.The Global Navigation Satellite System (GNSS) monitoring stations are primarily deployed along the section from Jinzhong Station to Taigu East Station, as illustrated in Fig. 6, which depicts the spatial distribution of 35 observation stations.Monitoring activities commenced in July 2021.

Regional deformation feature analysis
In the experiment, the SAR image captured on November 29, 2021, was selected as the primary reference image.Ground Control Points (GCP) generated by PS-InSAR were utilized for orbit refinement and re-flattening, thereby improving the accuracy of SBAS-InSAR processing 22 .Figure 5 illustrates the deformation rate results obtained by processing the research area using SBAS-InSAR and overlaying them on Google Earth Pro (V7.3.6.9345;software available at https:// www.google.com/ earth/ about/ versi ons/# downl oad-pro).The histogram's peak in the upper right corner indicates that a majority of the area experiences no deformation, which aligns with the actual conditions of the selected area, providing evidence for the reliability of the SBAS-InSAR monitoring results utilizing GCP points generated by PS-InSAR.The most significant subsidence area is located at the intersection of the two railways, exhibiting an average subsidence rate of -38 mm/year.In contrast, the main uplifted areas are primarily concentrated in the southwest region, specifically the Taigucheng District, and the villages in the northwest region, particularly in the south of Liujiabao Township.These uplifted areas have average subsidence rates of 15 mm/year and 20 mm/year, respectively.
There are three distinct deformations observed along the studied railway area, as depicted in Fig. 6.These subsidence areas are marked as A, B, and C, while three InSAR subsidence feature points are selected in each area and labeled as 1, 2, and 3, respectively.The selection of InSAR settlement feature point 1 in each area is determined using the nearest neighbor method with respect to the GNSS monitoring station near the railway.The time series of settlement measurements for each feature point in every region are presented in Fig. 7.The time series analysis reveals that the settlement trends and cumulative subsidence of the three feature points within area A, located at different positions, are primarily consistent.This suggests a large subsidence range within area A, corroborating the findings depicted in Fig. 5. Notably, characteristic point A-1, situated near the center of the Figure 8a represents the time series data after the coordinate conversion of the GNSS data using Eq. ( 16).It showcases the unified coordinate frame for the results monitored by the GNSS base stations near feature points A-1, B-1, and C-1, as well as the InSAR results.The results show that the Mean Absolute Error (MAE) of InSAR monitoring settlement data of feature point A-1 is 4.61 mm, at feature point B-1 is 3.83 mm, and at feature point C-1 is 2.28 mm.To enhance the persuasiveness of the results, we selected the remaining GNSS monitoring points and employed the nearest neighbor method to choose the corresponding InSAR points to validate the accuracy of the InSAR results.Due to the constraints of our research scope and our specific focus on areas of interest, there were situations where InSAR results were unavailable near some GNSS stations or the distance to the nearest InSAR feature point exceeded a certain threshold, rendering a direct comparison unfeasible.As a result, we selected the locations and monitoring data of 16 eligible GNSS points and compared them with nearby InSAR feature points using the nearest neighbor method.In Fig. 8b, the time series results of the coordinate transformation are presented for 16 randomly selected GNSS monitoring points along the railway and the adjacent InSAR feature points.
Table 2 shows the MAE and Root Mean Square Error (RMSE) of InSAR points near different GNSS stations.We can find that except for point 3, the MAE and RMSE of the InSAR points are both small.These results affirm the reliability and precision of SBAS-InSAR technology for monitoring settlement data.Therefore, this paper considers the InSAR monitoring data along the railway as the real observation data and conducts retrospective prediction research.

Experiment Experimental scheme
In this study, we utilized three prediction models: the ARMA model, the Prophet model, and the EEMD-Prophet prediction method, to predict and compare the InSAR monitoring points along the high-speed railway and its surrounding areas.Firstly, we utilized GNSS monitoring points along the railway in this area and applied the nearest neighbor method to select the corresponding InSAR settlement point.Subsequently, we conducted a comparative analysis of the settlement trend and the accuracy of InSAR monitoring.
Secondly, we utilized the three prediction models to forecast and compare the feature points along the railway.Initially, we predicted and analyzed the feature points located close to the selected GNSS monitoring points from the previous step and compared the results with the corresponding GNSS monitoring data.Subsequently, for the feature points situated at a greater distance from the GNSS monitoring station, we employed a retrospective prediction method for comparative analysis.
Thirdly, at a regional scale, we selected all InSAR points near Taigu East Railway Station along the railway line.Using the retrospective prediction method, we employed the three prediction methods to generate land subsidence predictions and created a subsidence map.We then compared this map with the original InSAR land subsidence map.
The BeiDou Navigation Satellite System (BDS) measurement and monitoring points deployed along the railway line are utilized to evaluate the accuracy of surface subsidence monitoring using through InSAR technology.BDS monitoring station can measure both horizontal and vertical deformations.The displacement of the BDS base station can be converted into the Line-of-Sight (LOS) direction of InSAR 23 .This conversion is represented by Eq. ( 15): In Eq. ( 15), d LOS represents the satellite line of sight direction, α h is the satellite radar heading angle, θ inc is the satellite radar incident angle, d N , d E , and d U are the north, east, and displacement in the upward direction.Assuming that the horizontal displacement of the BDS station is negligible 24 , it can be simplified as follows: To assess the prediction accuracy of feature points, we employ RMSE and MAE as evaluation metrics.The two parameters are defined as follows: where Ẑ is the predicted value, Z is the measured value.RMSE ∈ (0, +∞] , the closer the matter is to 0, the more consistent the expected value is with the measured value.MAE ∈ (0, +∞] , the closer the value is Closer to 0, then it reflects the combination model is better.www.nature.com/scientificreports/Among them, Z t x i , y i represents the predicted settlement value at time t and position x i , y i , while Z I,t x i , y i represents the observed settlement value obtained from InSAR monitoring at time t and position x i , y i .The coefficient of determination, R 2 , quantifies the degree to which the independent variable can explain the changes in the dependent variable.A value closer to 1 indicates a stronger fit of the predictive model.ẐI,t x i , y i represents the true value of the observed data, this article refers to the InSAR observed value.Z I,t x i , y i is the mean of the observed data.Z t is the predicted value of the model.

Single point forecasting and analysis
The subsidence time series data of feature points obtained by SBAS-InSAR technology are used as actual values for retrospective prediction research.Among them, the monitoring data from the initial 35 periods were selected as the sample data to predict the settlement for the next five periods, and the actual monitoring results were compared with the predicted results.
Firstly, EEMD decomposition is performed on the time series data of SBAS-InSAR monitoring points, resulting in the extraction of multiple IMF components and a residual component from the time series signal.This process allows us to capture the local characteristics of the monitoring points at different times, enabling the interpolation and prediction of their time series.By subjecting InSAR time series data to EEMD, the original signal can be approximated, and its components, referred to as IMF, can be extracted.Figure 9 presents a schematic diagram of the EEMD-decomposed IMF for a specific timing monitoring point of SBAS-InSAR.The time  www.nature.com/scientificreports/series acquired through SBAS-InSAR technology is decomposed into five components, with each IMF component reflecting distinct information and containing oscillation characteristics of the original sequence.IMF1 represents the component with the highest frequency, while the frequency decreases sequentially for the subsequent IMF components.IMF5, being the component with the lowest frequency, captures the trend information of the original sequence.The downward trend observed in IMF5 aligns with the overall trend of the original sequence.Table 1 clearly shows the SAR satellite has encountered data loss during a specific time interval, resulting in gaps in the subsidence time series used for SBAS-InSAR analysis of surface subsidence.These missing images present significant challenges to accurately assessing subsidence.In this study, the selected subsidence period spans from August 8, 2022, to October 7, 2022.Leveraging the periodicity of Sentinel-1A remote sensing images acquired every 12 days, we can infer the missing image timestamps by referencing to the information presented in Table 1.Specifically, the temporal gaps correspond to August 20, 2022, and September 25, 2022.Based on this premise, the initial selection for verification of the EEMD-Prophet prediction method includes feature points A-1, B-1, and C-1, which are in proximity to the GNSS base station.These selected feature points serve as the basis for conducting prediction experiments using the proposed EEMD-Prophet method.The predicted values for the missing time intervals are then compared with the corresponding GNSS measurements during those periods.
Figure 10 showcases the settlement prediction of feature points A-1, B-1, and C-1 utilizing the EEMD-Prophet method.It is evident that the projected values generated by the EEMD-Prophet method closely align with the GNSS measured values during the initial three prediction periods.However, the accuracy of the predictions gradually diminishes as the forecasting horizon extends.By incorporating span interpolation for the prediction time, the EEMD-Prophet method effectively anticipates the missing image data of SBAS-InSAR.The forecast results and corresponding relative error statistics for feature points A-1, B-1, and C-1 through the implementation of the EEMD-Prophet method are presented in Table 3. Notably, the missing images in InSAR correspond to August 20, 2022, and September 25, 2022.
To validate the effectiveness of the proposed method for settlement prediction along the railway, a single prediction method using EEMD-Prophet is insufficient.In this study, we conducted prediction experiments using the Prophet model, the ARMA model 5,6 , and the EEMD-Prophet model, as proposed in this paper.The experimental results were compared and analyzed to assess the feasibility of the proposed method.
For the purpose of comparative analysis, we selected feature points A-2, A-3, B-2, B-3, C-2, and C-3.These points were chosen because they are situated at a significant distance from the GNSS base station, and the accuracy of SBAS-InSAR technology meets the required standards.Consequently, the InSAR monitoring values of these feature points were employed as the actual observation data, and a retrospective prediction approach was adopted for comparative analysis and verification.The last five periods were designated as the forecast timeframe and compared against the actual observed values.
The prediction results depicted in Fig. 11

Regional forecast and analysis
In Sect.4.2, we employed the EEMD-Prophet, Prophet, and ARMA models to predict the settlement time series at a single feature point.In this section, we focus on the prediction of settlement at Taigu East Station, a significant area along the Taijiao high-speed railway.The geographical location of Taigu East Station is depicted in Fig. 4 Furthermore, Fig. 12 displays the region encompassing Taigu East Station and its 455 corresponding InSAR monitoring points in Google Earth Pro (V7.3.6.9345;https:// www.google.com/ earth/ about/ versi ons/# downl oad-pro).These monitoring points are utilized for prediction and evaluation using the EEMD-Prophet, Prophet, and ARMA models.Figure 13 is a visual representation of ground settlement generated through Kriging interpolation 26 after predicting the monitoring points in the Taigu East Station area.Upon examination, two prominent uplift areas and one notable subsidence area are identified within the vicinity of Taigu East Station.The subsidence area is situated at the intersection of the railway line and Taigu East Station.This observation suggests that, once trains commence operations after railway construction is completed, they exert pressure on the foundation, resulting in ground subsidence.Of particular note, the proximity of the railway station to the railway line indicates that the train load on this site may be substantial, resulting in noticeable settlement.
Figure 13 presents the settlement prediction of the Taigu East Station area using different prediction models, resulting in a ground settlement map.The EEMD-Prophet prediction model exhibits consistency with the actual InSAR land subsidence.However, as the prediction time increases, the established subsidence map gradually deviates from the actual InSAR ground subsidence map.To assess the accuracy and precision of each prediction method for regional prediction, we utilize the RMSE zone , coefficient of determination R 2 , and the time required for Taigu East Station prediction.Table 5 displays the average indicators of the prediction models at three prediction time points.
When examining the results, it becomes evident that the EEMD-Prophet prediction model outperforms both the ARMA model and the Prophet model in terms of regional prediction accuracy, achieving an RMSE zone of 3.

Conclusions
This paper utilizes 40 SAR images to monitor surface subsidence in the study area through the SBAS method.By analyzing the 40-period time-series subsidence changes and cumulative subsidence, the SBAS monitoring values are employed as the time-series data sample for retrospective forecasting research.For forecasting research, the EEMD-Prophet forecasting model, the Prophet model, and the ARMA model are employed.Overall, the combined model demonstrates higher prediction accuracy compared to the Prophet model and the ARMA model.
In the study of SBAS-InSAR settlement, this paper introduces the EEMD-Prophet model, which combines EEMD decomposition and Prophet prediction, to interpolate and predict InSAR monitoring time series.The proposed method exhibits strong adaptability and provides favorable prediction results for analyzing similar subsidence time series signal data.However, it should be noted that the decomposition and prediction process of the EEMD-Prophet model introduces an additional workload.Throughout the research process, it is crucial to consider that settlement changes along the high-speed railway are influenced by a variety of factors.Therefore, there is a need for further investigation into the development of an efficient and self-adaptive prediction method capable of capturing sudden changes over time.
In addition to the methods explored in this study, it's worth noting that machine learning techniques, such as LSTM, have shown promise in the prediction of ground subsidence.LSTM, with its ability to capture complex temporal dependencies in data, can be particularly useful in modeling and predicting land subsidence time series, especially in urban areas and regions influenced by anthropic activities.While the EEMD-Prophet model demonstrated enhanced predictive accuracy, combining it with machine learning methods could potentially offer even more robust results, especially in scenarios where sudden changes or non-linear patterns in subsidence need to be captured in a timely manner.This integration of techniques could represent an exciting avenue for further research in the field of land subsidence analysis and prediction.

( 2 )
Perform EEMD decomposition on the nonlinear time series results to obtain IMF components, and use the Prophet model to predict each IMF component separately.(3) Superimpose the prediction results of each IMF with equal weight to obtain the prediction results of the settlement point.(4) Evaluate and analyze the prediction accuracy of the settlement point.

−
k+α(t) T δ t− l+α(t) T γ close to the GNSS station, exhibits a maximum cumulative settlement of 110 mm.Feature points A-2 and A-3, situated at the periphery of the subsidence zone, display cumulative subsidence of 90 mm and 80 mm, respectively.During the monitoring period of the three feature points in area A, the cumulative settlement trend demonstrates a sharp increase, followed by a gentler rate, and then another sharp increase.The period of gentler cumulative settlement for the three feature points spans from August 2021 to February 2022, indicating a consistent subsidence trend along the railway lines in area A. The railway section in area B encompasses both uplifted and subsided regions.Among these, feature point B-1 experiences significant cumulative subsidence of 75 mm, while feature point B-2, situated between B-1 and B-3, demonstrates cumulative subsidence of 25 mm.Feature point B-3 represents an uplifted area with a relatively minor cumulative uplift.Area C is predominantly characterized by uplift, albeit at a slower rate compared to the sinking trend observed in area A. The maximum cumulative uplift within this area is observed at feature point C-1, reaching 22 mm, while featuring points C-2 and C-3 exhibit cumulative uplifts of 10 mm each.

Figure 5 .
Figure 5. Deformation rate and histogram of the study area (the color red indicates subsidence, the color blue indicates an uplift and The color green represents relatively stable regions).

Figure 6 .
Figure 6.Distribution of GNSS stations along the railway and areas A, B and C.

Figure 7 .
Figure 7. Settlement of feature points in ABC region.
23 mm.Compared to the ARMA model and the Prophet model, the EEMD-Prophet model yields a improve in RMSE zone by 43.53% and 12.47%, respectively.Moreover, the coefficient of determination for the EEMD-Prophet model is closer to 1, reaching 0.79, surpassing the other two prediction models.While the EEMD-Prophet model requires an additional 22s for EEMD decomposition on top of the Prophet model, its running speed exhibits a 12.45% improvement compared to the calculation time of the ARMA model.

Figure 12 .
Figure 12.Distribution of InSAR monitoring points at Taigu East Station.

Figure 13 .
Figure 13.Land subsidence of Taigu East Station regional prediction model.

Table 1 .
The acquisition time of SAR image data from 2021 to 2022.

Table 2 .
InSAR accuracy evaluation near GNSS points at different locations.

Table 3 .
Predictive value and relative error of EEMD-Prophet combined model.

Table 4 .
Prediction accuracy of different models.

RMSE/mm MAE/mm RMSE/mm MAE/mm RMSE/mm MAE/mm
The average RMSE and MAE values of the EEMD-Prophet method are recorded as 2.41 mm and 1.47 mm, respectively.The RMSE improves by 58.01% and 32.3% compared to the ARMA model and the Prophet model, respectively.Similarly, the MAE exhibits improvements of 62.69% and 33.78% when compared to the ARMA model and the Prophet model, respectively.

Table 5 .
Evaluation indicators of the model in regional prediction.